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Abstract. We measure the shear viscosity in a two-component Fermi gas of atoms, 
tuned to a broad s-wave coUisional (Feshbach) resonance. At resonance, the atoms 
strongly interact and exhibit universal behavior, where the equilibrium thermodynamic 
properties and the transport coefficients are universal functions of the density n and 
temperature T. We present a new calibration of the temperature as a function of global 
energy, which is directly measured from the cloud profiles. Using the calibration, 
the trap-averaged shear viscosity in units oi hn is determined as a function of the 
reduced temperature at the trap center, from nearly the ground state to the unitary 
two-body regime. Low temperature data is obtained from the damping rate of the 
radial breathing mode, while high temperature data is obtained from hydrodynamic 
expansion measurements. We also show that the best fit to the high temperature 
expansion data is obtained for a vanishing bulk viscosity. The measured trap-averaged 
entropy per particle and shear viscosity are used to estimate the ratio of the shear 
viscosity to the entropy density, which is compared that conjectured for a perfect fluid. 
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1. Introduction 

The measurement of the shear viscosity is currently of particular interest in the context 
of a recent conjecture, derived using string theory methods, which defines a perfect 
normal fluid [I]. The perfect fluid conjecture states that the ratio of the shear viscosity 
rj to the entropy density s has a universal minimum, 

1 > ^ A. (1) 

s - AnkB ^ ^ 

One example of a nearly perfect fluid is the quark-gluon plasma produced in a collision 
between two gold ions at an energy of 100 GeV per nucleon, which is thought to be a good 
approximation to the state of matter that existed microseconds after the Big Bang [2J. 
The collision produces a cigar-shaped plasma at a temperature of 2 x 10^^ K that exhibits 
elliptic flow, where the narrow direction of the plasma expands faster than the long 
direction, as observed in the momentum distribution [3]. A second example is a strongly 
interacting two-component Fermi gas of atoms, at a temperature of 10^^ K. Released 
from a cigar-shaped optical trap, a strongly interacting Fermi gas also exhibits elliptic 
flow, which is directly observed in the spatial profile of the expanding cloud Despite 
a difference in temperature of 19 orders of magnitude and a difference in a density of 25 
orders of magnitude, both systems exhibit nearly frictionless hydrodynamics and have 
a similar ratio r]/s. 

Ultracold strongly interacting Fermi gases are generally of broad interest, as they 
provide a tunable tabletop paradigm for strongly interacting systems, ranging from 
high temperature superconductors to nuclear matter. First observed in 2002, quantum 
degenerate, strongly interacting Fermi gases are being widely studied [11 [5l El [7]. To 
obtain strong interactions, characterized by a divergent s-wave scattering length [8], a 
bias magnetic field is used to tune the gas to a broad collisional (Feshbach) resonance, 
where the range of the collision potential is small compared to the interparticle spacing. 
In this so-called unitary regime, the two-body interaction potential produces no relevant 
length scales. Hence, the local thermodynamic and transport properties of a resonantly 
interacting gas are determined by the interparticle spacing L and the thermal de Broglie 
wavelength At, i.e., they are universal functions of the density n and temperature T [9]. 
Unitary Fermi gases therefore provide a scale-invariant system for exploring minimum 
viscosity hydrodynamics. 

The rj/s ratio is experimentally accessible in a trapped universal Fermi gas, 
where the entropy and other thermodynamic properties have been measured both 
globally [ini HU [121 1131 El] most recently, locally [iHl HE]. As the viscosity can 
be determined from hydrodynamic experiments [171 [THl [191 EO] , the predicted minimum 
ratio can be directly compared to that from Fermi gas experiments [T9l [201 121] . 

The scale of the rj/s ratio for a unitary Fermi gas can be understood using 
dimensional analysis. As noted above, there are only two natural length scales 
/ G {L, At}- Shear viscosity has units of momentum/area. As the natural momentum 
is of order h/l and the natural area is rj oc h/l^. At temperatures well below the 
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Fermi temperature at which degeneracy occurs, the Fermi momentum sets the scale so 
/ ~ L, and rj cc h/L^ cc hn. For a normal fluid above the critical temperature, the scale 
of entropy density s c:inkB and rj/s ~ h/ks- For much higher temperatures above the 
Fermi temperature, one expects that / ~ oc /iT~^/^, so that the shear viscosity scales 
as r/ oc h/X'^T oc T^'^/Ti^- 

2. Universal Hydrodynamics 

To properly measure the shear viscosity with high precision over a wide temperature 
range, we employ two experimental methods [21]. To determine the shear viscosity at 
low temperatures, we measure the damping rate of the radial breathing mode excited 
by releasing the cloud for a short time and then recapturing it [22| . For measurement at 
high temperatures, we observe the aspect ratio of the cloud as a function of time after 
the cloud is released from a deep optical trap. To consistently extract the viscosity from 
these two experiments, we use universal hydrodynamic equations, which contain both 
the friction force and the viscous heating rate. 

We determine the shear viscosity rj in the normal fluid regime by using a 
hydrodynamic description of a single component fluid, where the velocity field v(x, t) is 
determined by the scalar pressure and the shear viscosity pressure tensor, 

m{dt + v ■V)vi = fi + }^ diUtrap, (2) 

j ^ 

where f = —VP/n is the force per particle arising from the scalar pressure P, Utrap 
is the potential of the optical trap, and m is the atom mass. The second term on the 
right describes the friction forces arising from both shear rj and bulk viscosities, where 
aij = dvi/dxj + dvj/dxi — 25jjV ■ v/3 and a[j = 6ijV ■ v. 

For a unitary gas, the evolution equation for the pressure takes a simple form, since 
P = 2S/3 [HI [23], where £ is the local energy density (sum of the kinetic and interaction 
energy). Then, energy conservation and Eq. [2] implies [21] 

+ V ■ V + 5V ■ v/3)P = 2g/3. (3) 

Here, the heating rate per unit volume q = rj '^ij/'^ + C(V ■ v)^ arises from friction 
forces due to relative motion (shear viscosity) and dilation (bulk viscosity) of the volume 
elements. 

The evolution equation for the force per particle fi is readily determined by 
differentiating Eq. |3]with respect to Xi and using the continuity equation for the density, 

(5* + V ■ V + • v) /, + Y.i^^v,)f, - ^ (9,V ■ v) ^ = (4) 

Force balance in the trapping potential f/jrap(x), just before release of the cloud, 
determines the initial condition /i(0) = diUtrapi'^)- 

These hydrodynamic equations include both the force and the heating arising from 
shear and bulk viscosities. The solution is greatly simplified when the cloud is released 
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from a deep, nearly harmonic trapping potential Utrap, as /i(0) is then linear in the 
spatial coordinate. If we neglect viscosity, the force per particle and hence the velocity 
field remain linear functions of the spatial coordinates as the cloud expands. Thus 
c}j(V ■ v) = and the pressure P does not appear in Eq. HI Numerical integration [2l] 
shows that non-linearities in the velocity field are very small even if the viscosity is not 
zero, because dissipative forces tend to restore a linear flow profile. Hence, the evolution 
equations [2] and H] are only weakly dependent on the precise initial spatial profile of P 
and independent of the detailed thermodynamic properties. 

For the expansion and breathing mode experiments, we assume that the velocity 
field is exactly linear in the spatial coordinates. We take the force per particle to be of 
the form /j = ai{t)xi and assume that the density changes by a scale transformation [25], 
n(x, t) = n[x/bx{t),y /by(t), z/bz(t)]/T. Here T = b^byb^ is the volume scale factor and 
/ (i^xn(x, t) = is the total number of atoms, which is conserved. Current conservation 
then requires that the velocity field take the form Vi = bi(t)/bi(t), so that the pressure 
term in Eq. |4] vanishes as discussed above. 

The bulk viscosity is predicted to vanish in the normal fluid phase at unitarity [261 
[27] . so we initially exclude the bulk viscosity in our universal hydrodynamic 
equations [2] and m to extract the shear viscosity. In general, the (static) shear viscosity 
takes the universal form 

?7(x, t) = a{9) hn, (5) 

where 9 = T/Tpin) is the local reduced temperature and Tpin) = (Sn'^n)'^/^ / {2mkB) 
is the local Fermi temperature. Kinetic theory shows that rj ^ in the low density 
region of the cloud [28], as required for energy conservation. 

With these assumptions, the evolution equation for the force per particle, Eq. H] 
yields 



^^bi ^2 b, \ 2 



where, for zero bulk viscosity (C = 0), the heating rate is q = rj J2ijO'ij/2. Here, aij is 
evaluated using djVi = Sijbi/bi, which is spatially constant. Multiplying both sides by 
Xi and integrating over all space, the left side yields / d^iLnxf = N{xf) = N{x^)Qbj(t), 
where {x'^)o is the equilibrium mean square size of the trapped cloud in the i*'* direction. 
Integrating by parts, and assuming r/ — )■ as n — )■ 0, the right side is proportional to 
/ (i^xr7(x, t), yielding 

Similarly, Eq. |2] for the velocity field takes the form 

k m mN{x^)obKt)''''' 
Here, the a;? term on the right arises from the harmonic trapping potential, which is 
retained for the breathing mode and set equal to zero for expansion, where the cloud is 
released from the trap. 
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The right hand sides of Eqs. and [7] depend on the trap-averaged viscosity 
parameter, a, where 

a = ^ f d^:>^r]{:>^,t). (8) 



Nh 

The spatial integral exists, since as discussed above, as — )• as the density goes to 
zero at the edges of the cloud [2B] • As the viscosity produces a first order perturbation 
to perfect hydrodynamics, we can evaluate 6 in Eq. [5] using a zeroth order adiabatic 
approximation, so that 6 has a zero convective derivative everywhere. Since the number 
of atoms in a volume element is conserved along a stream tube, a is then time- 
independent. 

For measurements of the viscosity at high temperatures, the cloud is released from a 



deep optical trap and the aspect ratio ax{t)/az{t) = ^ {x^) / {z^) is observed as a function 
of time after release. Fig. [H Here, x and z refer to the initially narrow and long directions 
of the trapped cloud, respectively. Eqs. [6] and [7] are solved using the measured trap 
frequencies and the initial mean square cloud sizes, with the initial conditions 6j(0) = 1, 
hiifS) = 0, and aj(0) = muf. The ratio cr^(t)/o"j;(t) = (Uz/ux) bx(t)/bz(t) is determined 
as a function of time and compared to the data, yielding very good fits. Fig. [H with a 
as the only free parameter, which is determined by minimizing x^- 

We measure the viscosity at low temperature using a breathing mode [22], which 
is excited by a brief release and subsequent recapture. The oscillation of the transverse 
radius of the trapped cloud is observed as a function of time after excitation and the 
damping rate is measured by fitting an exponentially damped sinusoid to the data. 
For the breathing mode, the amplitude is small and the cloud radii change by a scale 
transformation of the form 6j = 1 + ti, with << 1. As the heating term in Eq. [6] 
containing afj is oc ef, the heating rate is negligible for the breathing mode and the 
force per particle evolves adiabatically to first order in small quantities, so that Eq. [6] 
yields = muf / (fe^F^/^). Using this in Eq. [71 we determine the breathing frequency ub 
and the damping rate 1/r. For a cylindrically symmetric trap, where Ux = ojy = uj±, 



ub = w± y 10/3. The damping rate arises from the an term, which, for a cylindrically 
symmetric trap is proportional to €x = iy = e±. 

More generally, we assume a nearly cylindrical trap potential as used in the 
experiments, with Ux, ojy the transverse trap frequencies and S = (ux—cOy) / ^uJxOJy « 1. 
Then, ub = ^JlOuxUJy/3. In terms of the cloud size observed for the x-direction, the 
damping rate is given by 

1 ^, . , 

The transverse mean square size of the trapped cloud before excitation is given by 
(a;^)o = (x^)/6^.(t), where (x^) is determined by imaging at a time t after the cloud is 
released and bx{t) is the calculated hydrodynamic expansion factor. We self-consistently 
determine bx{t) and a from the measured damping rate and the measured cloud size after 
expansion, using Eq. |9]as a constraint. Since the viscosity is small at low temperatures, 
we initially calculate bx{t) for zero viscosity to obtain an initial guess for (x^)o. This 
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yields an initial approximation for a. This initial value is then used in Eqs. El and [7] (with 
uf = and including the heating rate) to determine a better approximation to b^{t), 
which in turn yields a better value for (x^)o and a. This procedure quickly converges. 

3. Experiments 

Our experiments employ a 50-50 mixture of the two lowest hyperfine states of ^Li 
fermions, which is confined in a stable CO2 laser trap [29]. The mixture exhibits 
a Feshbach resonance at 834 G for which the s-wave scattering length diverges |30j . 
producing a unitary gas where the two-body scattering cross section is inversely 
proportional to the relative kinetic energy. 

After forced evaporation by lowering the trap depth, the trap is recompressed. At 
the final trap depth, parametric resonance is used to measure the oscillation frequencies 
of weakly interacting atoms. The frequencies obtained from the measurements are 
corrected for anharmonicity to determine the harmonic oscillation frequencies for 
energies small compared to the trap depth. For the deep optical trap used for the high 
temperature measurements, we obtain = 2tt x (182.7 ±0.5) Hz, = 2tt x (5283 ±10) 
Hz, ooy = 2tt X (5052 ± 10) Hz, and Co = {uj,^ujyLo,y/^ = 27r x (1696 ± 9) Hz. The total 
number of atoms ranges from = 4.0 x 10^ at E = 2.3 Ep to N = 6.0 x 10^ at 
E = 4.6 Ep, where less evaporation is employed. For A^ = 6.0 x 10^, the Fermi energy 
of an ideal gas at the trap center is Ep = (3A^)^/'^ hco = ks x 9.9 /iK, small compared to 
the trap depth Uq = ks x 460 /iK. 

For the high temperature experiments, the aspect ratio ax(t)/az(t) is measured 
as a function of time after release to determine a for energies E between 2.3 Ep and 
4.6 Ep, Fig. [H We also take expansion data at one low energy point E = 0.6 Ep, where 
the viscosity is small compared to that obtained at higher temperatures and the density 
profile is approximately a zero temperature Thomas- Fermi distribution. The black curve 
shows the fit for zero viscosity and no free parameters. To obtain a high signal to 
background ratio, we measure the aspect ratio only up to 1.4. For comparison, the 
green dashed curve shows the prediction for a ballistic gas. In these experiments, where 
the total energy of the gas E larger than 2 Ep, we find that the density profile is well fit 
by a Gaussian n(a;, 2;, t) = no{t) exp(—a;^/cr^ — ?/^/cr^ — z^/cr^), where (Tj(t) = 6j(t)o"j(0) 
is a time dependent width, no(t) = N / {-K^/'^axO'yaz) is the central density, and A^ is the 
total number of atoms. 

For the breathing mode experiments, the gas is initially cooled to nearly the 
ground state, where the total number of atoms is A^ = 2.0 x 10^, and then the 
trap is recompressed. Energy is then added by release and recapture, the gas is 
allowed to equilibrate, and the breathing mode is excited. The damping rate is 
measured for energies E between 0.5 E'j? and 2.5 Ep. For the shallow trap used in 
the breathing mode experiments, the measured parametric resonance frequencies are: 
ool = ^/ooJJ^ = 27T X 1696(10) Hz, ujuy = 1.107(0.004), and = 27t x 71(3) Hz, 
so that u) = {cOxCOyCOzy^^ = 27r x 589(5) Hz is the mean oscillation frequency and 
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A = u)z/u)± = 0.045 is the anisotropy parameter. The typical Fermi temperature 
Tp = (3NY^'^hu}/kB of a corresponding noninteracting gas is ~ 2.4 /iK, and the trap 
depth is Uo/kB = 35 /xK. The transverse mean square size (x^) is measured after an 
expansion time t = 1 ms. 

(A) 




T 1 1 1 1 1 1 1 1 

200 400 600 800 1000 1200 1400 1600 



Time After Release (us) 

(B) 

Figure 1. Anisotropic expansion. (A) Cloud absorption images for 0.2, .3, 0.6, 
0.9, 1.2 ms expansion time, E — 2.3 Ep', (B) Aspect ratio versus time. The expansion 
rate decreases at higher energy as the viscosity increases. Sohd curves: Hydrodynamic 
theory with the viscosity as the fit parameter. Error bars denote statistical fluctuations 
in the aspect ratio. 

The initial energy per particle E is measured from the density profile of the trapped 
cloud along any one direction by exploiting the virial theorem [HI |23] , which holds in 
the unitary regime. For measurements of the axial density profile, 

E = 3mul{z')[l-K{z')/al,], (10) 

where k = 15Ef/{4Uo), corrects for anharmonicity in the trapping potential |14j . 
Here {z"^) is the mean square size of the axial density profile of the trapped cloud 
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Figure 2. Trap-averaged viscosity coefficient a — J <P:s.r]/ {fiN) versus initial energy 
per atom. Blue circles: Breathing mode measurements; Red squares: Anisotropic 
expansion measurements. Bars denote statistical error arising from the uncertainty in 
E and the cloud dimensions. 



and apz = y'^Ep / {mul) is the Fermi radius for the z-direction. The energy also can 
be measured from the mean square size in the transverse direction, i.e., z — )■ x. A 
systematic uncertainty of 3% in Ep arises from the < 10% uncertainty in the absolute 
atom number N [31j . 

As shown in Fig. [H the expansion data are very well fit over the range of energies 
studied, using a as the only free parameter. We find that the friction force produces 
a curvature that matches the aspect ratio versus time data. For the expansion data, 
the indirect effect of heating is significant in increasing the outward force. Including 
heating in the analysis significantly increases the fitted a compared to that obtained 
when heating is omitted |21], as discussed below. For the breathing mode experiments 
at low temperature, the effect of heating is small, as discussed above. 

Together, the breathing mode and expansion measurements determine the fitted 
viscosity coefficients a for the entire energy range, as shown in Fig. [2l As shown in § [5l 
by calibrating the temperature, we determine a as function of reduced temperature at 
the trap center, prior to release. These results for the trap-averaged viscosity coefficient 
as a function of reduced temperature can be used to test predictions [32l |33l [31] , within 
the local density approximation. 

Despite the large values of a at the higher energies, the viscosity causes only a 
moderate perturbation to the adiabatic expansion, as shown by the expansion data 
and the fits in Fig. [H The breathing mode data and expansion data smoothly join, 
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Figure 3. Effect of tlie lieating rate in Eq. 2] on the measured viscosity 
coefficient a versus initial energy per atom. Solid(open) circles/squares: heating is 
included(excluded). Blue circles: low energy breathing mode data; Red squares: high 
energy expansion data. The high and low temperature data smoothly join only when 
heating is included. Bars denote errors arising from the uncertainty in E and the cloud 
dimensions. 



provided that the heating rate is included in the analysis. In contrast, omitting the 
heating rate produces a significant discontinuity, nearly a factor of 2, between the high 
and low temperature viscosity data, as shown in Fig. [31 The agreement between these 
very different measurements when heating is included shows that hydrodynamics in the 
universal regime is well described by Eqs. [2] and HI 

4. Vanishing of the Bulk Viscosity at Unitarity 

The bulk viscosity is predicted to vanish in the normal fluid phase at unitarity [261 l2Z] , 
which is why we did not include the bulk viscosity in our initial analysis to extract the 
shear viscosity. Energy conservation also requires the bulk viscosity to vanish in the 
two-body collision limit. We show that by including the bulk viscosity in both heating 
and force equations, the best fit to our anisotropic expansion data is the one for which 
the bulk viscosity exactly vanishes. 

By including the bulk viscosity in Eqs. El and HI with the same forms as used above 
for the density profile, the velocity field and the force per particle, we can apply a fit 
to our anisotropic expansion data. As shown in Fig. HI not only is the reduced of the 
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Figure 4. Comparison between tlie best fit with pure shear viscosity and the best 
fit with pure bulk viscosity. The green curve is the no free parameter prediction 
with zero shear and bulk viscosity. The black dots are the anisotropic expansion data 
at E/Ep — 3.3. Error bars denote statistical fluctuations in the aspect ratio. The red 
curve is the best fit with pure shear viscosity coeflicient as = 24.4 and reduced 
— 1.6. The blue curve is the best fit with pure bulk viscosity coefficient as = 16.7 
and reduced ~ 8.6. 



pure bulk viscosity fit much larger than that of the pure shear viscosity fit, but also it 
does not produce the curvature that our experimental data shows. This curvature arises 
from the fact that the anisotropic shear viscosity pressure tensor slows the transverse 
expansion and speeds up the axial expansion. In contrast, the bulk viscosity pressure 
tensor symmetrically slows the expansion in all directions. 

We also choose the bulk viscosity and apply a fit lo our data using the shear 
viscosity as the only free parameter [35j . As shown in Fig. |5l zero bulk viscosity gives the 
minimum reduced x^, which is consistent with the prediction of vanishing bulk viscosity 
in the normal fluid phase at unitarity [26l 



5. Shear Viscosity versus Reduced Temperature 

For comparison with predictions for the temperature dependence of the shear 
viscosity [HI [321 [331 [M] , we give the trap-averaged viscosity coefficient a as a function 
of reduced temperature at the trap center, prior to release of the cloud. 
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Figure 5. P^r degree of freedom versus bulk viscosity with shear viscosity as the 
only free parameter. 



Here, the local Fermi temperature at the trap center, before the cloud is released, 
is given by Tp(?T,o) = /i^(37r^no)^''^/(2mfcs) and Tpi = Ep/kB = Tpini) is the ideal 
gas Fermi temperature at the trap center, with nj the ideal gas central density for a 
zero temperature Thomas- Fermi distribution, SN/^Tr'^aFxO'FyO'Fz)- As shown previously, 
{nj/riQy/^ oc E/Ep. Hence, at higher temperatures, where Tq/Tp oc E/Ep, we have 
^0 oc (E/Epf [21]. 

To measure Tq/Tpj for the high temperature expansion data, where E > 2Ep, we 
use the second virial coefficient approximation to the local energy density |36] . 



S = {'i/2)nkBT{l + B2n\%) 



(12) 



where n is the total density and At = h/y/2jTmkBT is the thermal wavelength. 
B2 = 1/2''/^— 62/"\/2 is the second virial coefficient for a unitary gas, with 62 = 1/2, which 
is universal and known to be accurate for experiments in this temperature regime [T31IT6] . 
Here, the first term in B2 arises from degeneracy for each spin state of the 50-50 
mixture and the second term arises from interactions between them. Force balance 
in the trap requires [231 E] / d^^S = J d^xx ■ VUtrap{^)/2. For the Gaussian density 
profile observed in the high temperature experiments, one obtains 



at 
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(13) 



Note that Eq. [13] includes first order corrections arising from B2 and trap anharmonicity. 
For the high temperature data, where the density profile is Gaussian, we have 
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Figure 6. Total energy per particle of a strongly interacting Fermi gas at 840 G versus 
the entropy per particle. The blue dots show the entropy obtained by adiabatically 
sweeping the magnetic field from 840G to 1200G and using an exact many-body 
calculation |13j for the entropy at 1200 G, where kpa — —0.75. [See 6*1200 Table 
1, Ref. [S]]. The red dots are the theoretical calculations using the second virial 
coefficient approximation including the trap anharmonicity for the trap condition of 
Ref. [llj . The green curve is the power law fit with a discontinuous heat capacity, as 
described in Eq. [151 

{nj/uQy^^ = 4:{al/apJ/7i^^^. Using the initial Tq/Tfj obtained from the cloud profile, 
we determine 6q from Eq. [11] 

We find Tq/Tfi for the low temperature breathing mode data from our previous 
measurements of energy versus entropy, for which the trap is as shallow as that used 
for the low temperature breathing mode measurements. The low temperature E versus 
S data employs the entropy S'J'goo given in Table 1 of Ref. [H], which is corrected for 
the finite interaction strength at 1200 G [37]. To properly calibrate the temperature 
with high precision over a wide range, we join the the experimental E{S) data with 
theoretical calculations of E and S in the high temperature regime by exploiting the 
second virial coefficient approximation for a unitary Fermi gas, as shown in Fig. [6] For 
simplicity, we fit E{S) data using a smooth curve with a discontinuous heat capacity [H] 
as follows: 

Ey{S) = El + cS"^; S>Sc. (14) 

Constraining the values of Ei and c by demanding that energy and temperature be 
continuous at the joining point Sc, one obtains: 

E^S) = Eo + aS''] 0<S<Sc 
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Figure 7. Measured energy versus the temperature obtained from tlie calibration of 
Eq. [15] (red dots) ; For comparison, we show the data obtained by the ENS group [T6] 
(black dots) and the theory of Hu et al. f38] (green curve). 

E>{S) = Eo+ aS^^[l - b/d + h/d{S/ScY]; S > S,. (15) 

We fix the ground state energy Eq = 0.47 by using (3 = —0.60 [H]. Then we use 
a fit, where we give equal weighting to the low temperature data points and the 
calculated high temperature points. We find a = 0.10(1), b = 1.57(15), d = 2.23(3). 
The critical parameters obtained from the fit are Sc = 2.04(39), E^ = 0.78(12) Ep. 
Then, T = dE/dS yields Tc = 0.24(8) Tpi- Note that the relatively large error bar for 
Sc and hence Tc arises from joining measured data with the calculated high temperature 
points, and fitting Ey{S) by a single power law. However, this method yields a smooth 
temperature calibration that reproduces the temperatures used in the virial calculation 
to better than 3%, Fig. [6l Then, we obtain E/Ep as a function of Tq/Tpi, which 
yields Tq/Tpj from the measured initial energy of the cloud. Fig. [71 The central density 
no, prior to release is determined from column density of the trapped cloud by fitting 
the spatial profile with a Gaussian distribution, which is adequate except at the lowest 
temperature where the density profile is a zero-temperature Thomas-Fermi Distribution. 

Fig. [8] shows the trap-averaged viscosity coefficient a versus the initial reduced 
temperature at the trap center 6q, from nearly the ground state to the unitary two- 
body regime. The inset in Fig. [8] is the predicted local shear viscosity in the normal 
fiuid regime versus the local reduced temperature [Mj . 
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Figure 8. Trap-averaged viscosity coefficient a — J d^xri/{hN) versus reduced 
temperature 9q ~ Tq/TpIuq) at the trap center, prior to release. Blue dots: Breathing 
mode measurements; Red dots: Anisotropic expansion measurements. Bars denote 
statistical errors arising from the uncertainty in a, E and the cloud dimensions. 
The green curve shows the fit a = ^3/2^0^^ + 0:1/2^0^^ with 0:3/2 = 2.96(3) and 
"1/2 = 1.87(8), for the temperature from nearly the superfluid transition point up 
to the two-body unitary regime. The inset shows the theoretical prediction for the 
local shear viscosity a in the normal fluid phase versus the local reduced temperature 
9 = T/Tpin) from Ref. [31]. The red curve shows the fit a = 03/2 6*3/2+01/2 6*1/2 ^-^^-^ 
03/2 = 2.11 and ai/2 = 0.74. 



5.1. Temperature Scaling 

In the universal regime, the local viscosity at high temperatures is expected to scale as 
T^/^, as described in § [1] The viscosity at the trap center prior to release of the cloud 
can be written as r/o = aohuQ. In the high temperature limit, we then have 

"0 = «3/2 6'o'^^ (16) 

where 03/2 is a universal coefficient pjj. As 9 has a zero convective derivative everywhere 
(in the zeroth order adiabatic approximation), ^0 at the trap center where v = 0, has 
a zero time derivative and is therefore constant as is the trap-averaged viscosity 
coefficient a. 
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We note that in the hydro dynamic regime, the local value of oc T3/2 in the hieh 
temperature limit is independent of density and spatially constant, so that a formally 
does not exist. To test the prediction of the T'^/^ temperature scaling in the high 
temperature regime, we assume that r] relaxes to the equilibrium value in the center of 
the trap, but vanishes in the low density region so that a is well defined. As noted 
in § [21 this behavior is predicted by kinetic theory [28] and is required for energy 
conservation. Assuming a ~ ao, we observe the predicted T^/^ temperature scaling 
in the high temperature regime [21]. In this case, fitting the data with a = 03/2^0^^ 
yields 03/2 = 3.4. This is consistent with predictions, where kinetic theory [39] shows 
that the trap-averaged T^/^ coefficient should be larger by a factor of ^ 1.3 than that 
of the local value for the two-body unitary regime, 0:3/2 = 2.77 [TBJ. However, more 
work is needed to fully understand the relation between the local viscosity and the 
trap-averaged viscosity parameter measured in the experiments. 

In Fig. [HI both the experimental data and the theoretical calculations are found 
to scale as T^/^ in the high temperature regime, consistent with Eq. [161 However, we 
observe that the damping rate of the radial breathing mode reaches a plateau at the 
higher temperatures [201 [22]. This flattening can be explained if a has a term oc E, 
since 1/r oc a/E, according to Eq. [9] Since 6^^^ oc E at higher temperatures, we fit the 
data of Fig. [8] with the two-parameter fit function, 

a = a3/2 6l^'^ + (17) 

where a is the trap averaged viscosity and 6*0 is the local reduced temperature at the 
trap center. Fitting all of the data for the normal fluid regime {E > 0.7 Ep), we find 
'S3/2 = 2.96(3) and ai/2 = 1.87(8), where the errors are statistical from the fit. In this 
case, 03/2 = 2.96 is closer to the local value predicted for the high temperature two-body 
unitary regime, where 03/2 = 2.77. 

We can also use the same two-parameter fit for the predicted local viscosity in the 
normal fiuid regime [3l] , inset Fig. [HI 

a = a3/2^'/' + ai/2^'/', (18) 

where a is the local shear viscosity and 9 = T/Tpin) is the local reduced temperature. 
We obtain 03/2 = 2.11 and ai/2 = 0.74. 

The coefficient 0:1/2 = 1.87 obtained from the data is much larger than the local 
value obtained from the prediction 01/2 = 0.74. However, to compare these parameters, 
we need the trap-average of the reduced temperature, (0^/^) = 1.84 6'^^. Hence, the fit 
to the calculated viscosity yields O1/2 = 1.84 01/2 = 1.84 x 0.74 = 1.36, somewhat lower 
than the fit to the data. However, if we scale the fit coefficients for the predicted viscosity 
by a factor 2.77/2.11 = 1.31, then 03/2 agrees by construction with the accepted high 
temperature two-body limit, while the predicted O1/2 increases to 1.31 x 1.36 = 1.78, in 
good agreement with the measurements. 

In the low temperature normal fiuid regime, Fermi liquid theory [IQ] predicts that 
the local viscosity should have an upturn oc 1/6"^. After trap- averaging, we find 

the predicted viscosity near the critical temperature is about a factor of 2 larger than 
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the observed value, with a temperature dependence that is not consistent with our data. 
In earher work |18] , pair formation was suggested as a means of avoiding Pauh blocking, 
which would otherwise cause the viscosity to increase with decreasing temperature. 
Recent predictions [32] suggest that the onset of a normal state pairing gap and the 
formation of pairs may suppress the viscosity at low temperature, which is consistent 
with our observations. 

5.2. Ratio of the shear viscosity to the entropy density 

To estimate the ratio ?7/s, we use the approximation 77/5 = ahn/s = {h/kB)oi/{s/nkB) — 
{h/kB)a/S, where 5" is the average entropy per particle of the trapped gas in units of 
ks- For the low temperature data, we obtain 5* from the fitting curve in Fig. El 

To determine S for the high temperature data, we employ the second virial 
coefficient approximation of Eq. [121 Assuming that the density has Gaussian profile, as 
observed in the experiments and used to find the energy, we obtain the trap-averaged 
entropy per particle in terms of the reduced temperature at the trap center. 



Then we calculate E versus Oq (using Eqs. [lOl [HI and[l3]) to obtain E versus S. Fig. M 
shows the average entropy per particle over the measured energy range. 

Fig. [10 shows the estimated ratio rj/s. The inset shows the low temperature 
behavior, which is about 5 times the string theory limit (red dashed line) near the 
critical energy [Hj Ec/Ep = 0.7 — 0.8. This result is in good agreement with recent 
predictions [H], where r]/s = 4.7. We note also that the apparent decrease of the rj/s 
ratio as the energy approaches the ground state 0.48 Ep [13| does not require that the 
local ratio — >■ as T — )■ 0, since contributions from the cloud edges significantly increase 
S compared to the local s at the center. 

6. Summary 

We have measured the trap-averaged shear viscosity coefficient as a function of reduced 
temperature at the trap center, from nearly the ground state to the two-body unitary 
regime. From our measurements of global entropy and energy, we have calibrated 
the temperature versus energy and find very good agreement with integrated local 
measurements as well as predictions. Using the temperature calibration, we determine 
the trap-averaged shear viscosity coefficient as a function of reduced temperature at the 
trap center, which is compared to recent predictions. We show that the best fit to our 
experimental expansion data is obtained with a vanishing bulk viscosity. The measured 
trap-averaged entropy per particle and shear viscosity are used to estimate the ratio of 
the shear viscosity to the entropy density. Near the transition point, the ratio is found 
to be ~ 5 times that of a perfect fiuid. More work is needed to understand the behavior 




(19) 
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Figure 9. Calculated trap- averaged entropy per particle versus average energy per 
particle for the trapped gas. Blue: Low temperature data calculated from the fit in 
Fig. |6l where error bars arise from the energy uncertainty. Red: High temperature 
calculation for the deep trap used in the expansion experiments, based on the second 
virial coefficient. Error bars arise from the energy uncertainty. 

of the viscosity in the low temperature regime, where the breathing mode damping rate 
and hence the viscosity appear to approach as T — )■ 0. 
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